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Abstract 

Measurements are presented of the CP violation observables S and C in the decays 
of and mesons to the J/'ipKg final state. The data sample corresponds to 
an integrated luminosity of 3.0 fb“^ collected with the LHCb experiment in proton- 
proton collisions at center-of-mass energies of 7 and 8 TeV. The analysis of the time 
evolution of 41 500 B^ and B^ decays yields S = 0.731 ± 0.035 (stat) ± 0.020 (syst) 
and C =—0.038 ± 0.032 (stat) ± 0.005 (syst). In the Standard Model, S equals 
sin(2/3) to a good level of precision. The values are consistent with the current world 
averages and with the Standard Model expectations. 
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The violation of charge-parity (CP) conservation in processes involving B mesons was 
hrst observed in the “golden mode” P° —)■ J/iIj by the BaBar and Belle experiments at 
the asymmetric e^e“ colliders PEP-II and KEKB [^|^. Since then, measurements of CP 
violation in this decay mode have reached a precision at the level of 10“^ |^[^. Thus, these 
measurements play an important role in constraining and testing the quark-flavor sector 
of the Standard Model [^|^, which relates CP-violating observables to a single irreducible 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [^|^. As the 
J/ijjKg hnal state is common to both the P° and the P° meson decays, the interference 
between the amplitudes for the direct decay and for the decay after oscillation 

results in a decay-time dependent CP asymmetry between the time-dependent decay rates 
of B^ and P° mesons. 




S sin( Am t) — C cos( Am t) 


r(P0(f)^ J/il’KI) + r(P0(f)^ cosh(^) + AArsinh(^) 




( 1 ) 


Here, B^{t) and B^{t) indicate the flavor of the B meson at production, while t indicates 
the decay time. The parameters Am and AP are the mass and the decay width differences 
between the heavy and light mass eigenstates of the P°-P° system, and S, C, and Aaf 
are CP observables. As AP is negligible for the P°-P° system j^, the time-dependent 
asymmetry simplihes to A{t) = S'sin(Amt) — C cos(Amt). 

The P° —>■ J/'ifjKg decay is dominated by a 6—ccs transitionj^ and CP violation 
in the decay is expected to be negligible at the current level of experimental precision, 
giving C ^ 0. This allows to identify S with sin(2/3), where /3 = 

is one of the angles of the CKM triangle. Other measurements that constrain this 
triangle predict sin(2/3) as 0.771±q^II 10 , giving a small discrepancy with respect to 


the average of direct measurements, 0.682 ± 0.019 j^, where the most precise input 
comes from a CP violation measurement in P°—)■ J/ipK^ decays by the Belle experiment, 
S = 0.670 ± 0.029 (stat) ± 0.013 (syst) |^. To clarify the CKM picture, both better 
experimental precision and improved understanding of higher-order contributions to the 


decay amplitudes are required 11,12 


The analysis presented in this Letter supersedes a previous measurement by LHCb 13 
which was performed on data corresponding to an integrated luminosity of 1.0 fb“^ at 
a center-of-mass energy of 7TeV. By adding data corresponding to 2fb“^ at 8TeV and 
using an optimized selection and additional “flavor tagging” algorithms to identify the 
quark content of the B meson at production, we increase the statistical power of the 
analysis by almost a factor 6. 


The LHCb detector 14,15 is a single-arm forward spectrometer covering the pseu¬ 


dorapidity range 2 < r/ < 5, designed for the study of particles containing 6 or c quarks. 
The detector includes a high-precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet, and three stations of silicon-strip detectors and straw drift 


^Mention of a particular decay mode implies the inclusion of charge-conjugate states except when the 
measurement of CP violation is involved. 
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tubes placed downstream of the magnet. Different types of charged hadrons are distin¬ 
guished using information from two ring-imaging Cherenkov detectors. Photon, electron, 
and hadron candidates are identihed by a calorimeter system consisting of scintillating-pad 
and preshower detectors, an electromagnetic calorimeter, and a hadronic calorimeter. 
Muons are identified by a system composed of alternating layers of iron and multiwire 
proportional chambers. The online event selection system (trigger) consists of a 
hardware stage, based on information from the calorimeter and muon systems, followed by 
a software stage. 

The analysis is performed with —)■ J/'ip Kg candidates reconstructed in the J/^jj —)■ 

li~ and Kg —)■ 7r’''7r“ hnal states. Two oppositely charged particles, identified as muons 
with high momentum and high transverse momentum, are required to originate from 
a common space-point (vertex) and to have an invariant mass in a range ±60 MeV/ 
around the known J/xf) mass Since the meson has a lifetime of 1.5 ps and has high 
momentum, the resulting J/'ij) candidate is required to be significantly separated from all 
reconstructed pp collision points (primary vertices (PVs)). The Kg candidates are formed 
from two oppositely charged, high-momentum pion candidates with a clear separation 
from any PV in the event. Candidates decaying early enough for the final-state pions 
to be reconstructed in the vertex detector are characterized as long candidates and are 
required to have an invariant mass within ±15 MeV/ of the known Kg mass 17 . The 


Kg candidates that decay later, such that track segments of the pions cannot be formed 
in the vertex detector, are called downstream candidates; these have a larger momentum 
resolution than the long candidates, and thus the corresponding 7r’''7r“ pairs are required 
to have an invariant mass within ±55MeV/c^ of the known Kg mass. A good vertex 
fit quality and sufficient separation from the B^ decay vertex are required for the Kg 
candidate’s decay vertex. To eliminate background contributions from 7l°—)■ J/'ipA decays, 
the 7r~^ (vr”) candidate has to fulfill particle identification requirements if the invariant 
mass under a pn~ (vr+p) mass hypothesis is compatible with the A mass. 

The B^ candidates are reconstructed from J/'il) and Kg candidates that form a good 
quality vertex. Their decay time t is obtained from a fit to the full decay chain while 
constraining the production vertex of the B^ to the associated PV 18 . The reconstructed 


B^ candidate mass m is obtained from a similar fit with the pApT and 7r’''7r“ invariant 
masses constrained to the known J/il) and Kg masses. The latter fit must satisfy loose 
requirements on its quality, and resulting candidates are only retained if 5230 < m < 
5330 MeV/and 0.3 < t < 18.3 ps. The fit uncertainty at on the decay time is required 
to be smaller than 200 fs, which is well above the average resolution of 55 fs (65 fs) for 
candidates with long (downstream) Kg daughters. The quantity at is used later in the 
analysis as an estimate of the per-candidate decay-time resolution. Multiple PVs and, in a 
small fraction of events, multiple 5°—>■ J/'il)Kg candidates, lead to multiple (i?°,PV) pairs 
per event. In such cases, one pair is chosen at random. 

Various simulated data samples are used in the analysis. In the simulation, pp collisions 
are generated using Pythia [T^ with a specific LHCb conhguration 20 . Decays of hadronic 


particles are described by EvtGen 21 . The interaction of the generated particles with 
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the detector, and its response, are implemented using the Geant 4 toolkit 
in Ref. 


23 


22 as described 


Tagging algorithms are used to infer the initial flavor of the B meson candidate, i.e. 
whether it contained a 6 or a 6 quark at production. Each algorithm provides a decision d 
on the flavor of the B meson candidate (tag), and an estimate rj of the probability for that 
decision to be incorrect (mistag probability). The knowledge of the B meson production 
flavor is essential for this analysis, and so only candidates for which the tagging algorithms 
yield a decision are considered. 

One class of flavor tagging algorithms, the opposite-side tagger (OS), exploits the 
dominant production mechanism of b hadrons, i.e. the production of bb quark pairs, by 
reconstructing the b hadron produced in association with the signal B meson. The OS 
tagger uses the charge of the electron or muon from semileptonic b decays, the charge of 
the kaon from the b^ s decay chain, and the inclusive charge of particles associated 
with the secondary vertex reconstructed from the b hadron decay products; further details 
are described in Ref. 24 . 

is the inclusion of the same-side 


A major improvement in this analysis over Ref. 13 


pion tagger (SStt), which deduces the production flavor by exploiting pions produced in the 
fragmentation of the b quark that produced the signal B meson or in the decay of excited 
B mesons into the signal B meson 25,26 . Tagging pion candidates are selected requiring 


charged, high momentum and high transverse-momentum particles that are consistent 
with originating from the associated PV. Pions are identihed using information from the 
particle identihcation detectors, and the difference between the invariant mass of the B 
and the Btt^ pair is required to be less than 1.2 GeV/c^. Additionally, the flight directions 
of the pion and the B candidate must be compatible. If multiple pion candidates pass the 
selection, the one with the highest transverse momentum is used. The mistag probability 
is obtained using a neural network which is trained on simulated events and whose inputs 
are global event properties and kinematic and geometric information on the pion and B 
signal candidates. 

The tagging calibration is performed in control samples of B mesons whose hnal state 
determines the B flavor at decay time, by determining a linear correction u:{rf) that relates 
the estimated mistag probability r] with the mistag probability u observed in the control 
sample. To account for asymmetries in the detection efficiency of charged particles, which 
can lead to different mistag probabilities for B^ and B^ mesons, an additional linear 
correction function Au{r]) is determined. Asymmetries in the efficiency of the algorithms 
in determining a decision are found to be negligible. 

The B^ —)■ J/'il) decay is used to determine the flavor tagging calibration for the 
OS tagger. A consistency check of the calibration is performed in a control sample of 
B^ —)■ J/'ip K*^ decays, showing a good correspondence of the calibration between B^ and 
B^ decays. As the quarks that accompany the b quark in B~^ and B^ mesons differ, the SSvr 
tagger calibration is performed with )■ J/'ipK*^ decays 27 . Systematic uncertainties 


are assigned for the uncertainties associated with the calibration method and for the 
validity of the calibration in the signal decay mode. A summary of the calibration results 
is given in [the Appendix 
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Figure 1: Distribution of (a) the reconstructed mass and (b) logarithmic distribution of the decay 
time of tagged J/ipK^ candidates. The solid black lines show the fit projections, while the 

dashed (dotted) lines show the projections for the signal (background) components only. 


The effective tagging efficiency is the prodnct of the probability for reaching a 
tagging decision, ^tag = (36.54 ± 0.14) %, and the sqnare of the effective dilntion, 
D = 1 — 2u} = (28.75 ± 0.24) %, which corresponds to an effective mistag probability of 
oj = (35.62 ± 0.12) %. Compared to the previous LHCb analysis 13 the effective tagging 
efficiency Ses = increases from 2.38 % to 3.02 %, mainly due to the inclusion of the 

SSvr tagger. 

The values of the CP violation observables S and C are estimated by maximizing the 
likelihood of a probability density function (PDF) describing the unbinned distributions of 
the following observables: the reconstructed mass m, the decay time t and its uncertainty 
estimate at, the OS and SStt flavor tag decisions dos and dssvr, and the corresponding per- 
candidate mistag probability estimates r^os and r^ssjr. The fit is performed simultaneously 
in 24 independent subsamples, chosen according to data-taking conditions (7TeV, 8TeV), 
Tfg type (downstream, long), flavor tagging algorithm (OS only, SSvr only, OS and SSvr), 
and two trigger requirements. In each category the data distribution is modeled using 
a sum of two individual PDFs, one for the signal and one for the combinatorial 
background. 

The reconstructed mass of the signal component is parametrized with a double-sided 
Hypatia PDF with tail parameters determined from simulation. An exponential 
function is used to model the background component, with independent parameters for the 
downstream and long subsamples. The ht to the mass distributions yields 41560 ± 270 
tagged )■ J/%IjK^ signal decays. The mass distribution and projections of the PDFs 
are shown in Fig. [^(a). 

The decay-time resolution is modeled by a sum of three Gaussian functions with 
common mean, but different widths, which are convolved with the PDFs describing the 
decay-time distributions. Two of the widths are given by the per-candidate resolution 
estimate at, each calibrated with independent linear calibration functions. The third 
Gaussian describes the resolution for candidates associated to a wrong PV. The scale and 
width parameters are obtained in a ht to the decay-time distribution of a control sample of 
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candidates formed from prompt J/'ip and mesons. The parameters are determined 
separately for candidates formed from downstream and long ifg candidates. 

Trigger, reconstruction, and selection criteria distort the measured decay-time dis¬ 
tribution, leading to a decay-time dependent efficiency. Effects of the trigger requirements, 
which distort the decay-time distribution at low decay times, are determined using data. 
The misreconstruction of tracks leads to inefficiencies at large decay times. To account for 
this effect, an additional decay-time dependent efficiency of the form is used, where 
Pt is obtained from simulation. 

The PDF of true decay times t' is given by 


{t' 1 dosi dssAvoSi Vsstt) 





(1 — d'Ap) e ^ {1 — d'S sin(Amt') -|- d'C cos(AmT)} , 


( 2 ) 


where the tag decision d takes the value +1 (—1) for a tagged B^ (B^) candidate and d' 
takes the value -1-1 (—1) for the B^ component of the signal distribution, r is the B^ 
meson lifetime, and 


C{dj,Vj,d') = l + dj (^1-2 


Uj{T]j) + d‘ 


A^iVj 


( 3 ) 


represents the calibration of the tagging response from the tagging algorithm j = 
{OS,SS7r}. The production asymmetry Ap = [cr(i?°) — cr(S°)]/[(T(5°) -|-(T(i?°)], where 
a denotes the production cross-section inside the LHCb acceptance, is obtained us¬ 
ing a measurement in 7TeV pp collisions 29 . Considering differences between the 


7 and 8TeV data-taking conditions, the production asymmetries are determined as 
^7TeV ^ _o.oi08 ± 0.0052 (stat) ± 0.0014 (syst) and ^ ^ 

0.0004 ± 0.0018 (syst) [^. The background decay-time distribution is parametrized by 
a sum of exponential functions, convolved with the resolution model used for the signal. 
This parametrization does not depend on the tag decision and mistag probability estimates. 
The number of required exponential functions varies across subsamples. The decay-time 
distribution and projections of the PDFs are shown in Fig. [^(b). The distributions of the 
per-candidate resolution estimate at and the per-candidate mistag probabilities, r^os and 
r^ssjT, are modeled by empirical functions. Independent parameterizations are chosen for 
the signal and background components. 

The likelihood is a function of 83 free parameters, including S and C, and 48 yield 
parameters for the signal and the background components in 24 individual subsamples. 
Eleven parameters are external inputs, including the production asymmetry, the flavor 
tagging calibration parameters, and the mass difference Am 17 . These are constrained in 


the ht within their statistical uncertainties and taking their correlations into account. The 
likelihood ht yields S = 0.729 ± 0.035 and C = —0.033 ± 0.032 with a correlation coeffi¬ 
cient of p{S,C) = 0.483. Fig. shows the decay-time dependent signal-yield asymmetry. 
An additional ht with hxed C = 0 yields S = 0.746 ± 0.030. Corrections of -1-0.002 for S 
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Figure 2: Time-dependent signal-yield asymmetry (Ai^o—-/V^go)/(A^jjo + -/Vso). Here, (A^^jo) is 
the number of ^ J/tf^Kg decays with a {B^) flavor tag. The data points are obtained with 
the sPlot technique , assigning signal weights to the events based on a fit to the reconstructed 
mass distribution. The solid curve is the projection of the signal PDF. 


and —0.005 for C are applied to account for CP violation in mixing and for the 

difference in the nuclear cross-sections in material between and states 31 . The 
correction is negligible for the result for S with (7 = 0. 

Various sources of systematic uncertainties on the CP observables are examined, 
in particular from mismodeling PDFs and from systematic uncertainties on the input 
parameters. In each study, a large set of pseudoexperiments is simulated using a PDF 
modihed such as to include the systematic effect of interest; the relevant distributions 
from these pseudoexperiments are then htted with the nominal PDF. Signihcant average 
deviations of the ht results from the input values are used as estimates of systematic 
uncertainties. The largest systematic uncertainty on S, ±0.018, accounts for possible tag 
asymmetries in the background; for C the largest uncertainty, ±0.0034, results from the 
systematic uncertainty on Am. Systematic uncertainties on the flavor tagging calibration 
account for the second largest systematic uncertainty on S, ±0.006, and on (7, ±0.0024. 
The third largest uncertainty on S, ±0.005, arises from assuming AF = 0 and is evaluated 
by generating pseudoexperiments with AF set to the value of its current uncertainty, 
0.007ps“^ 1^, and then neglecting it in the ht. Remaining uncertainties due to neglecting 
correlations between the reconstructed mass and decay time of the candidates, mismodeling 
of the decay-time resolution and efficiency, the systematic uncertainty of the production 
asymmetry, and the uncertainty on the length scale of the vertex detector are small and are 
given in the Appendix Adding all contributions in quadrature results in total systematic 
uncertainties of ±0.020 on S and ±0.005 on C. 

Several consistency checks are performed by splitting the data set according to differ¬ 
ent data-taking conditions, tagging algorithms, and different reconstruction and trigger 
requirements. All results show good agreement with the nominal results. 

In conclusion, a measurement of CP violation in the interference between the direct 
decay and the decay after oscillation to a J/'ll) Kg hnal state is performed using 

41500 havor-tagged B^ —)■ J/ipKg decays reconstructed with the LHCb detector in a 


6 













sample of proton-proton collisions at center-of-mass energies of 7 and 8 TeV, corresponding 
to an integrated luminosity of 3.0 fb“^. The CP observables S and C, which allow the 
determination of the CKM angle /3, are measured to be 

S = 0.731 ± 0.035 (stat) ± 0.020 (syst), 

C = -0.038 ± 0.032 (stat) ± 0.005 (syst), 

with a statistical correlation coefficient p{S, C) = 0.483. When C is hxed to zero the 
measurement yields S = sin(2/3) = 0.746 ± 0.030 (stat). This measurement supersedes 
the previous LHCb result obtained with 1.0 fb“^ [^, and represents the most precise 
time-dependent CP violation measurement at a hadron collider to date. Furthermore, the 
result has a similar precision to, and is in good agreement with, previous measurements 
performed at the Belle and BaBar experiments at the KEKB and PEP-II colliders [^|^. 
This result is in excellent agreement with the expectations from other measurements and 
improves the consistency of the CKM sector of the Standard Model. 
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Appendix 


Overview of tagging calibration parameters 

The calibration functions of the mistag probability uj{t]) and the mistag probability 
difference Au{ri) = are chosen as 


=Pi{v- iv)) + Po, Au{r]) = Api{r] - {p)) + Apo . (Al) 

The OS calibration parameters are determined to be 

pOS ^ Q_3g^5 ^ Q_QQ33 ^ Q_QQ3g ('gyg^) ^ 

pOS ^ Q gyg ^ Q_Q 32 (stat) ± 0.009 (syst) , 

Ap°® = 0.0148 ± 0.0016 (stat) ± 0.0008 (syst) , (A2) 

ApOS = 0.070 ± 0.018 (stat) ± 0.004 (syst) , 

(ryOS^ = 0.3786 . 


Fig. ^ shows the raw mixing asymmetry and the calibration of the mistag estimates. The 
measured SStt calibration parameters are 


pf^ = 0.4232 ± 0.0029 (stat) ± 0.0028 (syst) , 

pf’^= 1.011 ±0.064 (stat) ±0.031 (syst), 

Apf^ = -0.0026 ± 0.0043 (stat) ± 0.0027 (syst) , (A3) 

Apf’^ =-0.171 ±0.096 (stat) ±0.04 (syst), 

(rySS^) = 0.425 . 


In the dataset of 114000 —)■ J/tl’Kg decays, the OS tagger algorithm yields a 

tagging power of (2.63 ± 0.04)% and the SSvr tagger algorithm yields (0.376 ± 0.024)%. 
For events that have both OS and SStt tagging decisions, the effective tagging power is 
(0.503 ± 0.010)%. The combined effective tagging power of these three overlapping tagging 
categories is (3.02 ± 0.05)%. 


Summary of systematic uncertainties 


The systematic uncertainties are summarised in Table Al 


certainty is calculated by summing the single uncertainties 
systematic uncertainties compared to the central values of S 
Here, we set S = 0.729 and C = —0.033 as reference. 


The overall systematic un¬ 
in quadrature. The relative 
and C are given in brackets. 
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Figure Al: (a) Raw mixing asymmetry (A^unmixed - Admixed)/(A^unmixed + Admixed) for all SStt 
tagged ;■ J/^pK*^ decay candidates. Here, A'unmixed (Admixed) is the number of ;■ J/pK*^ 
decays with a hnal state that does (not) correspond to the flavor tag. The black line shows the fit 
projection, (b) The linear calibration of the SSvr mistag probability with )• J/pK*^ decays. 


Table Al: Systematic uncertainties as and ac on S and C. Entries marked with a dash represent 
studies where no significant effect is observed. 


Origin 

(7s 


(7c 

Background tagging asymmetry 

0.0179 

(2.5%) 

0.0015 

( 4.5%) 

Tagging calibration 

0.0062 

(0.9%) 

0.0024 

( 7.2%) 

AT 

0.0047 

(0.6%) 


— 

Fraction of wrong PV component 

0.0021 

(0.3%) 

0.0011 

( 3.3%) 

z-scale 

0.0012 

(0.2%) 

0.0023 

( 7.0%) 

Am 



0.0034 

(10.3%) 

Upper decay time acceptance 



0.0012 

( 3.6%) 

Correlation between mass and decay time 




— 

Decay time resolution calibration 




— 

Decay time resolution offset 




— 

Low decay time acceptance 




— 

Production asymmetry 




— 

Sum 

0.020 

(2.7%) 

0.005 

(15.2%) 
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